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Sporadic  E  layers  in  the  lower  thermosphere  play  a  very  important  role  in  radio 
communications.  The  lifetimes  of  these  transient  layers  are  largely  controlled  by  the 
ion-molecule  chemistry  of  metallic  ions,  particularly  iron  (Fe),  that  are  produced  by 
meteoric  ablation.  The  aim  of  this  project  was  to  study  the  role  of  atomic  oxygen  in 
controlling  the  rate  at  which  Fe+  ions  are  neutralised.  The  rate  coefficients  for  the 
reactions  of  Fe+  with  N2O,  O2  and  N2,  and  the  reactions  of  FeO+  with  CO  and  O,  were 
measured  using  a  fast  flow  tube  coupled  to  a  mass  spectrometer,  with  Fe+  ions 
produced  by  the  pulsed  laser  ablation  of  an  iron  rod.  Very  good  agreement  was 
achieved  with  previous  measurements  of  the  first  four  reactions.  The  reaction  of 
FeO+  and  O  was  found  to  be  about  100  times  slower  than  the  Langevin  limit.  The 
results  were  then  employed  in  a  mesospheric  model,  which  is  now  able  to  predict  the 
top-side  scale  height  of  the  atomic  Fe  layer  and  the  ratio  of  Fe+/FeO+  above  90  km  in 
good  accord  with  observations.  Finally,  a  combination  of  emission  spectroscopy  and 
ab  initio  quantum  theory  was  used  to  show  that  FeO+  does  not  have  an  excited 
electronic  state  that  would  provide  a  suitable  LIDAR  probe  in  the  atmospheric 
transmission  window. 
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A.  Introduction 


Sporadic  E  layers  are  thin  layers  of  concentrated  plasma  that  occur  in  the 
lower  thermosphere  between  about  90  and  120  km.  They  play  a  very  important  role  in 
facilitating  over-the-horizon  radio  transmission  as  well  as  obstructing  space-to-ground 
communication,  depending  on  the  time-of-day  and  transmission  frequency.1  The 
objective  of  this  proposal  was  to  study  the  critical  ion-molecule  chemistry  that 
determines  the  lifetimes  of  these  transient  events. 

Measurements  made  with  rocket-borne  mass  spectrometers  have  shown  that 
sporadic  E  layers  are  largely  composed  of  the  metallic  ions  Fe+  and  Mg+,  produced  by 
meteoric  ablation.2  These  ions  (and  an  equal  concentration  of  electrons)  can  be 
concentrated  into  layers  only  a  few  kilometres  thick  by  various  mechanisms, 
including  wind  shear,  auroral  precipitation,  and  gravity  waves  within  the  polar  cap.3,4 
Above  about  1 10  km,  dielectric  recombination  (i.e.  Fe+  +  e"  — »  Fe  +  hv)  is  the  only 
significant  route  by  which  metallic  ions  can  be  removed.  However,  this  is  a  very 
inefficient  process  with  a  rate  coefficient  of  about  10" 1 1  cm3  s’1,5  corresponding  to  a 
lifetime  for  Fe+  of  about  1  day  even  in  a  very  intense  sporadic  E  layer  (JES  =10  MHz). 

In  contrast,  ion-molecule  reactions  become  important  below  110  km, 
especially  the  reaction 

Fe+  +  03  ->  FeO+  +  02  (1) 

We  have  recently  measured  the  rate  coefficient,  k\,  and  find  that  the  reaction  is  very 
rapid  below  100  km.6  The  resulting  molecular  ion  can  then  undergo  dissociative 
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recombination  with  an  electron: 

FeO+  +  e"  — >  Fe  +  O  (2) 

If  this  combination  of  reactions  1-2  were  able  to  occur  unchecked,  then  a  sporadic  E 
layer  descending  below  100  km  would  have  a  lifetime  of  only  a  few  minutes,  contrary 
to  observation.  In  practice,  atomic  O  must  play  a  crucial  role  through  the  reaction 

FeO+  +  O  -»  Fe+  +  02  (3) 

By  competing  with  reaction  2  this  limits  the  efficiency  with  which  Fe+  is  converted  to 
Fe,  and  hence  governs  the  lifetime  of  sporadic  E  in  the  90  -  110  km  region.  In  spite 
of  its  importance,  the  rate  coefficient  k 3  has  not  been  measured. 

Plane  et  al.7  have  estimated  fc3  from  detailed  balancing  with  the  reverse  reaction, 

Q 

whose  excitation  function  was  measured  in  a  molecular  beam  study.  This  showed 
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that  fc3  could  be  as  low  as  2  x  10"  cm  molecule"  s"  at  lower  thermospheric 
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temperatures  (200  K),  whereas  atmospheric  modelling  of  the  neutral  mesospheric  Fe 
layer  indicates  a  value  of  about  7  x  10" 12  cm3  s"1.  Simultaneous  measurement  of  Fe+ 
and  FeO+  in  the  lower  thermosphere2  also  indicate  a  value  for  /c2  around  10"11  cm3  s"1. 
It  is  clear  that  /c3  needs  to  be  measured  before  the  temporal  behaviour  of  sporadic  E 
can  be  properly  modelled;  this  was  the  major  focus  of  this  project. 

A  secondary  goal  was  to  study  the  laser  induced  fluorescence  (LIF)  spectrum  of 
FeO+  in  the  520  -  550  nm  region,  where  the  molecule  could  have  a  strong  6X  <—  X6Z 
transition,  by  analogy  with  the  isoelectronic  molecule  MnO.  A  recent  photofragment 
study9  has  shown  that  FeO+  and  MnO  are  strikingly  similar  molecules. 


B.  Description  of  the  Work 

Laboratory  Measurements 

Figure  1  illustrates  the  apparatus  employed  for  the  project.  A  stainless  steel 
fast  flow  tube  was  coupled  to  a  two-stage  quadrupole  mass  spectrometer  configured 
for  measuring  positive  ions.  The  source  of  Fe+  ions  was  the  pulsed  laser  ablation  of  a 
rotating  iron  rod  in  the  upstream  section  of  the  flow  tube,  using  a  frequency-doubled 
Nd:YAG  laser  (wavelength  =  532  nm,  energy  ~  20  mJ  pulse"1).  This  produced  pulses 
of  Fe+  (and  also  neutral  Fe)  that  were  entrained  in  a  substantial  flow  of  an  inert  gas 
(He),  and  swept  down  the  tube  at  velocities  in  excess  of  10  m  s’1.  The  pulses  of  ions 
were  then  detected  by  the  mass  spectrometer  operating  in  pulse-counting  mode,  and 
interfaced  to  a  multi-channel  scalar.  Figure  2  illustrates  a  sequence  of  Fe+  pulses 
recorded  at  different  flow  times,  demonstrating  the  loss  of  Fe+  by  diffusion  to  the 
walls  of  the  tube.  Because  the  pulse  shapes  are  near  Gaussian,  the  relative  Fe+ 
concentration  could  be  computed  either  from  the  pulse  maximum  or  the  integrated 
area  of  each  pulse.  The  advantage  of  working  with  pulses  is  that  the  flow  velocity  of 
the  ions  down  the  tube  can  be  measured  precisely,  whereas  when  using  a  continuous 
flow  of  ions,  as  in  most  conventional  flow  tubes,  the  plug  velocity  has  to  be  corrected 
to  account  for  the  wall  losses  which  create  strong  radial  gradients. 

The  reactants  (N2O,  O2,  O  etc.)  were  then  admitted  to  the  tube  using  either 
fixed  inlets  along  the  tube  or  a  sliding  injector.  The  latter  was  used  to  vary  the 
distance  (and  hence  time)  during  which  the  Fe+  and  a  reactant  were  in  contact. 

Atomic  O  was  generated  by  the  microwave  discharge  of  N2,  followed  by  titration  with 
NO  (N  +  NO  — >  N2  +  O).  The  absolute  concentration  of  O  in  the  flow  tube  was 
measured  by  the  classical  titration  with  NO2.  The  rate  of  loss  of  O  on  the  flow  tube 
walls  was  determined  by  adding  NO  just  upstream  of  the  optical  port  (see  Figure  1), 
and  observing  chemiluminescence  around  600  nm  as  the  flow  time  was  varied. 

In  order  to  search  for  possible  electronic  transitions  of  FeO+  in  the  visible  (500 
-  600  nm),  ozone  was  added  to  the  pulses  of  Fe+/Fe  just  before  the  optical  port.  The 
resulting  chemiluminescence  was  focused  into  a  fibre  optic  and  dispersed  onto  a  CCD 
detector  through  a  0.5  m  spectrometer  with  a  choice  of  three  gratings,  providing 
resolution  from  0.3  to  1.2  nm.  Because  strong  chemiluminescence  is  generated  from 
neutral  Fe  reacting  with  O3  (the  “orange  arc  bands”),  spectra  were  also  recorded  with 
N2O  added  upstream.  This  reacts  rapidly  with  Fe+  but  not  with  Fe  at  room 
temperature,10  providing  a  means  of  discriminating  between  emission  due  to  Fe+  and 
Fe. 

Experimental  Results 

In  order  to  characterise  this  novel  system  for  studying  ion-molecule  reactions, 
the  reactions  of  Fe+  with  N2O,  O2  and  N2  were  studied  and  compared  with 
measurements  made  previously  by  the  UEA  group  using  a  completely  different 
technique.6,10  As  shown  in  Table  1,  very  satisfactory  agreement  was  achieved. 

Before  studying  reaction  3,  we  also  decided  to  study  the  reaction 

FeO+  +  CO  -»  Fe+  +  C02 


(4) 


This  is  analogous  to  the  O  atom  reaction  3,  but  without  the  added  complication  of 
significant  wall  losses  of  the  reactant  down  the  tube.  FeO+  was  produced  by  the 
reaction 

Fe+  +  N20  -»  FeO+  +  N20  (5) 

Figure  3  illustrates  the  results:  Fe+  and  FeO+  were  monitored  by  the  mass 
spectrometer,  and  their  ratio  is  plotted  against  the  ratio  of  [N2O]  to  [CO].  This 
enabled  £4  to  be  obtained  directly  relative  to  the  rate  coefficient  for  Fe+  +  N2O,  which 
is  well  known  (Table  1  and  ref.  10). 

Reaction  3  was  then  studied  by  using  reaction  5  to  produce  FeO+  and  adding  O 
through  the  sliding  injector.  The  fraction  of  Fe+  relative  to  total  Fe+  (that  is,  in  the 
absence  of  either  N2O  or  O),  was  monitored  as  a  function  of  both  reactants.  Figure  4 
shows  an  experiment  where  the  N2O  concentration  was  varied  while  the  O  atom 
concentration  was  fixed.  This  illustrates  that  there  is  a  significantly  higher  fraction  of 
Fe+  in  the  presence  of  atomic  O,  confirming  that  reaction  3  does  indeed  occur.  In 
order  to  obtain  £3,  a  full  kinetic  model  of  the  flow  tube  was  developed,  which 
included  the  wall  losses  of  Fe+,  FeO+  and  O.  Four  data-sets  under  different  conditions 
were  then  optimised  to  yield  the  value  for  k 3  listed  in  Table  1.  An  example  of  the 
model  fit  is  shown  in  Figure  4. 

Figure  5  illustrates  a  chemiluminescence  spectrum  (blue  line)  produced  by 
adding  O3  to  the  Fe+/Fe  pulse.  When  N2O  was  added  to  remove  Fe+,  the  resulting 
spectrum  (red  line)  was  found  to  be  identical;  this  is  confirmed  by  the  absence  of 
structure  in  the  residual  spectrum  (dark  green  line,  values  multiplied  by  2  x  104). 
Indeed,  for  a  series  of  spectra  taken  over  the  range  400  to  900  nm,  the  observed 
chemiluminescence  was  shown  to  arise  exclusively  from  excited  states  of  FeO. 

Theoretical  Calculations 


The  experimental  search  for  electronic  transitions  of  FeO+  was  complemented 
by  a  set  of  ab  initio  quantum  calculations.  The  ground  state  of  the  ion  was  first 
optimised  using  the  B3LYP  hybrid  density  functional/Hartree  Fock  method  with  the 
6-311+g(2d,p)  basis  set  from  within  the  Gaussian  98  suite  of  programs.12  These 
reasonably  high  level  calculations  produced  excellent  agreement  with  measured 
parameters  (shown  in  parenthesis,  from  ref.  9):  bond  length  re  =  1.641  [1.648]  A, 
vibrational  frequency  ?  e  =  830  [838±4]  cm1,  bond  energy  Do(Fe-0+)  =  337  [335±5] 
kJ  mol1.  Calculations  were  then  performed  on  the  first  10  excited  states  of  FeO+, 
using  the  configuration  interaction  (Cl-Singles)  method.  “  The  vertical  excitation 
energies  and  transition  dipole  moments  are  listed  in  Table  2. 

In  order  to  check  the  ability  of  the  Cl-Singles  method  to  predict  accurate 
excitation  energies  for  these  transition  metal  oxides,  we  performed  analogous 
calculations  at  the  same  level  of  theory  on  MnO  (isoelectronic  with  FeO+)  and  FeO, 
since  the  spectroscopies  of  both  species  are  known.  The  results  for  these  molecules 
are  also  listed  in  Table  2. 

Discussion 


The  rate  coefficient  for  reaction  3  is  about  two  orders  of  magnitude  less  than 
the  typical  Langevin  limit  for  an  ion-molecule  reaction.  This  is  in  accord  with  the 
small  reaction  cross  section  for  the  reverse  reaction,  Fe+  +  O2  — >  FeO+  +  O,  which 
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peaks  at  only  1.7  A".  Indeed,  detailed  balancing  indicates  that  reaction  3  should 
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proceed  at  only  3x10"  cm  molecule"  s"  at  296  K,  a  factor  of  20  less  than  our 
experimental  measurement  (Table  1).  This  must  indicate  that  reaction  3,  which  is 
exothermic  by  158  ±  5  kJ  mol"1,  is  able  to  access  additional  channels  to  yield  products 
in  excited  states,  including  Fe+(4F  and  4D)  and/or  CM '  ?  g).  Although  we  have  so  far 
only  determined  k 3  at  room  temperature;  the  excitation  function  for  the  reverse 
reaction  indicates  that  reaction  3  should  have  a  very  small  temperature  dependence. 

Inspection  of  Table  2  shows  that  the  Cl-Singles  method  predicts  excited  states 
for  MnO  and  FeO  at  470  nm  and  579  nm,  respectively.  These  are  in  very  good 
accord  with  the  6Z  X6E  transition  of  MnO  at  520  nm,  and  the  5 A  X3A  transition 
of  FeO  around  580  nm  (Figure  5).  By  contrast,  the  first  excited  state  of  FeO+  with  a 
strongly  allowed  vertical  transition  is  4.7  eV  above  the  ground  state,  corresponding  to 
a  transition  in  the  UV.  This  result  is  supported  by  the  absence  in  the  spectrum  in 
Figure  5  of  chemiluminescence  produced  by  reaction  1,  which  is  sufficiently 
exothermic  to  populate  electronic  states  of  FeO+  up  to  2.6  eV  above  the  ground  state. 
With  no  evidence  for  a  suitable  transition  for  FeO+  in  the  visible,  it  was  decided  not  to 
attempt  a  laser  induced  fluorescence  scan. 

Atmospheric  Modeling 

Figure  6  compares  the  annual  average  Fe  layer  at  40°N  predicted  by  our  1- 
dimensional  model7  with  a  comprehensive  set  of  lidar  measurements  from  Prof.  C.  S. 
Gardner’s  group  at  U.  of  Illinois,  Urbana-Champaign.  The  scale  height  on  the 
topside  of  the  Fe  layer  is  a  sensitive  test  of  the  ion-molecule  chemistry  under 
consideration  here.  Excellent  agreement  was  achieved  by  using  the  measured  kj,  from 
this  study  and  fitting  £2,  the  rate  coefficient  for  the  dissociative  recombination  of  FeO+ 
with  an  electron.  The  best  fit  value  of  M  =  5  x  10"  cm  molecule"  s"  is  in  sensible 
accord  with  the  measured  rate  coefficients  for  the  analogous  reactions  of  NO+  and 
°2+- 

There  is  a  limited  set  of  simultaneous  observations  of  Fe+  and  FeO+  made  by 
rocket-bome  mass  spectrometry  in  the  lower  thermosphere.”  Prof.  E.  Kopp  (U.  of 
Bern)  has  made  this  rocket  data  available  to  us.  Figure  7  illustrates  profiles  of  the 
Fe+/FeO+  ratio  from  four  flights  where  FeO+  was  above  the  detection  limit.  The  ratio 
peaks  at  about  92  km,  and  is  between  150  and  250,  although  with  quite  large 
uncertainty  because  the  small  concentrations  of  FeO+  are  close  to  the  detection  limit. 
Our  model  currently  predicts  a  ratio  around  100  at  92  km.  This  is  in  sensible  accord 
with  the  rocket  data,  since  the  ratio  is  very  sensitive  to  the  O3/O  ratio  and  these 
species  have  significant  mixing  ratio  gradients  in  this  region  of  the  atmosphere. 

Finally,  the  dramatic  effect  of  this  chemistry  on  the  lifetime  of  Fe+  in  a 
sporadic  E  layer  is  illustrated  in  Figure  8.  This  shows  the  lifetime  of  Fe+  against 
being  neutralised  to  Fe,  as  a  function  of  height  and  electron  density.  For  critical 
frequencies  (JES)  ranging  from  1  to  10  MHz,  the  corresponding  electron  densities  in  a 
sporadic  E  layer  vary  from  about  1  x  104  to  1  x  106  cm"3.  Figure  8  shows  that  very 
intense  layers  (fEs  >  8  MHz,  [e  ]  >  8  x  105  cm"3)  will  have  lifetimes  of  about  1  hour  at 
100  km,  and  should  be  transient  below  this  altitude.  In  contrast,  weak  layers  (fEs  <  2 
MHz,  [e‘]  >  5  x  104  cm"3),  should  have  lifetimes  of  several  hours  down  to  92  km. 

Personnel 

The  personnel  involved  in  the  project  were  Dr.  Tomas  Vondrak  (senior  research 
associate),  Mr.  Kenneth  Woodcock  (graduate  student),  and  Prof.  John  Plane. 


Table  1.  Rate  coefficients  (296  K)  determined  during  this  project. 


Reaction 

Rate  coefficient  a 

This  study  Previous  work 

Fe+  +  N20  FeO+ 

(4.99  ±  0.24)  x  10'11 

(4.47  ±  0.94)  x  10‘llb 

Fe+  +  N2  +  He  -»  FeN2  +  He 

(2.37  ±  0.20)  x  10'30 

(2.10  ±  0.32)  x  10'30c 

Fe+  +  02  +  He  — »  FeO  2  +  He 

(3.54  ±  0.24)  x  10'30 

(4.17  ±  0.23)  x  10'30c 

FeO+  +  CO  — >  Fe+  +  C02 

(1.49  ±0.  06)  x  10'10 

(2.05  ±  0.2)  x  10'10d 

FeO+  +  0  Fe+  +  02 

(6  ±  2)x  10'12 

Not  measured 

a  Units:  bimolecular  reactions,  cm3  molecule1  s'1;  termolecular,  cm6  molecule  2  s'1 
b  Ref.  10.  c  Ref.  6.  dRef.  11 


Table  2.  Excited  states  of  FeO+,  MnO  and  FeO,  calculated  using  the  Cl-Singles 
method  with  the  6-311+g(2d,p)  basis  set.  The  vertical  excitation  energy, 
corresponding  wavelength  and  transition  dipole  moment  (f)  are  shown.  The 
highlighted  transitions  indicate  the  first  excited  state  of  each  molecule  with  a  strongly- 
allowed  transition  from  the  ground  state. 


FeO+  excited  states 

MnO  excited  states 

FeO  excited  states 

1 

0.21  eV  6048  nm  f=0.0000 

0.21  eV  5794  nm  f=0.0000 

0.052  eV  23979  nm  f=0.0000 

2 

0.30  eV  4185  nm  f=0.0000 

0.92  eV  1346  nm  f=0.0015 

0.052  eV  23979  nm  f=0.0000 

3 

0.38  eV  3289  nm  f=0.0000 

2.64  eV  470  nm  f=0.2234 

0.30  eV  4185  nm  f=0.0000 

4 

0.88  eV  1414  nm  f=0.0001 

2.75  eV  450  nm  f=0.2257 

1.00  eV  1238  nm  f=0.0000 

5 

1.09  eV  1136  nm  f=0.0000 

3.39  eV  365  nm  f=0.1786 

1.00  eV  1238  nm  f=0.0000 

6 

1.48  eV  837  nm  f=0.0000 

4.04  eV  306  nm  f=0.0010 

1.24  eV  1004  nm  f=0.0000 

7 

4.70  eV  263  nm  f=0.0056 

4.47  eV  277  nm  f=0.0000 

1.61  eV  769  nm  f=0.0001 

8 

5.20  eV  238  nm  f=0.0041 

4.47  eV  277  nm  f=0.0000 

1.61  eV  769  nm  f=0.0001 

9 

5.33  eV  232  nm  f=0.0012 

4.78  eV  260  nm  f=0.0233 

1.88  eV  661  nm  f=0.0000 

10 

5.37  eV  231  nm  f=0.0144 

5.59  eV  222  nm  f=0.0295 

2.14  eV  579  nm  f=0.0559 
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Pump 


Figure  1.  Schematic  diagram  of  the  fast  flow  tube/pulsed  laser  ablation  system  for 
studying  the  kinetics  of  ion-molecule  reactions  of  Fe-containing  species. 


Figure  2.  Time  profiles  of  Fe+  pulses  for  different  bath  gas  (He)  flows.  From  left  to 
right:  5500  seem,  5000  seem,  4500  seem,  4000  seem,  3500  seem  and  3000  seem. 
The  total  drift  length  is  89.5  cm. 
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Figure  3.  Study  of  the  reaction  between  FeO+  and  CO.  The  FeO+/Fe+  ratio  detected 
at  the  downstream  end  of  the  flow  tube  is  plotted  against  the  ratio  of  FTO  and  CO. 
The  slope  yields  the  rate  coefficient  relative  to  that  for  Fe+  +  N2O.  N2O  and  CO  were 
added  8.5  cm  and  13.5  cm  downstream  of  the  Fe+  source,  respectively.  The  total  drift 
distance  was  89.5  cm.  The  flow  tube  pressure  was  1.2  Torr  and  the  total  gas  flow  was 
3500  seem. 
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Figure  4.  Study  of  the  reaction  between  FeO+  and  atomic  O.  The  ratio  of  Fe+  in  the 
presence  of  N2O,  versus  the  total  Fe+  in  the  absence  of  reactants,  is  plotted  against  the 

1  o 

N2O  concentration  (black  points).  The  experiment  is  then  repeated  with  3.7  x  10 
cm" 3  of  atomic  O  added  (blue  points).  The  model  simulations  are  the  green  and  red 
points,  respectively. 


Wavelength  (nm) 

Figure  5.  Emission  spectra  recorded  when  O3  (3  x  1013  cm  3)  is  added  to  the  pulse  of 
Fe+/Fe  just  before  the  optical  detection  port  at  the  downstream  end  of  the  flow  tube. 
The  spectra  were  recorded  with  1.2  nm  resolution  on  a  CCD  detector:  (a) 
chemiluminescence  from  Fe  and  Fe+  +  O3;  (b)  chemiluminescence  from  Fe  +  O3,  with 
the  Fe+  removed  by  addition  of  N2O  upstream;  (c)  the  residual  (a)  -  (b),  multiplied  by 
2  x  104,  demonstrating  no  identifiable  chemiluminescence  due  to  Fe+  +  O3  producing 
excited  states  of  FeO+. 


Figure  6.  The  annual  average  layer  of  atomic  Fe  at  40°N,  comparing  lidar 
observations  from  the  University  of  Illinois  (green  line)  with  the  UEA  1-D  model 
prediction  (black  line).  The  predicted  profile  of  Fe+  (red  hne)  is  also  shown. 
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Figure  7.  The  ratio  of  Fe+/FeO+  measured  by  rocket-borne  mass  spectrometry  (Prof. 
E.  Kopp,  U.  of  Bern).  For  the  sake  of  clarity,  each  flight  is  shifted  by  100  units.  The 
rocket  flights  were  made  during  twilight  at  high  latitudes  during  summer. 
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Figure  8.  The  predicted  lifetime  of  Fe+,  with  respect  to  being  neutralised  to  Fe,  as  a 
function  of  height  and  electron  concentration.  Conditions:  Janaury,  midnight,  40°N. 


